Electron-hole asymmetry in two-terminal graphene devices 
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A theoretical model is proposed to describe asymmetric gate- voltage dependence of conductance 
and noise in two-terminal ballistic graphene devices. The model is analyzed independently within 
the self-consistent Hartree and Thomas-Fermi approximations. Our results justify the prominent 
role of metal contacts in recent experiments with suspended graphene flakes. The contact-induced 
electrostatic potentials in graphene demonstrate a power-law decay with the exponent varying from 
— 1 to —0.5. Within our model we explain electron-hole asymmetry and strong Fabri-Perot oscilla- 
tions of the conductance and noise at positive doping, which were observed in many experiments 
with submicrometer samples. Limitations of the Thomas-Fermi approximation in a vicinity of the 
Dirac point are discussed. 

PACS numbers: 72.80.Vp,73.23.Ad,73.40.Cg 
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I. INTRODUCTION 

The recent observation of the fractional quantum 
Hall effect in graphene 1 calls for better understanding 
of two-terminal transport measurements in suspended 
graphene devices. The two-terminal geometry has also 
been used to reach the ballistic regime of transport with 
submicrometer flakes on Si0 2 substrate. 2,3 These and 
other experiments 4 ' 5 share such common features as the 
electron-hole asymmetry in the gate- voltage dependence 
of the conductance and the prominent conductance os- 
cillations at positive doping that are beyond a simplistic 
theory of ballistic transport. 6 ' 7 In this paper we attribute 
these experimental observations to the effect of charge 
transfer in a vicinity of metal electrodes and to a weak 
screening in graphene at low doping. 

The influence of metal contacts on electron transport 
in nearly ballistic graphene has been investigated both 
experimentally 8-12 and theoretically. 13-19 An agreement 
between theory and experiment has, however, been lim- 
ited due to both the complexity of experimental setups 
and a lack of important ingredients in the considered 
theoretical models. In this study we suggest a mini- 
mal model, which leads to a very good description of 
experimental data by Xu Du et al. 5 on charge transport 
in suspended graphene. The gate-voltage dependence of 
transport properties predicted by the model also agree 
qualitatively with those observed by Hcerschc et al. 2 for 
graphene on Si0 2 substrate. 

The electrostatic potential landscape in graphene with 
metal contacts has been imaged by means of scanning 
photocurrent microscopy 8 ' 11 to demonstrate the long- 
range decay of the contact-induced potential (about 1 
/mi). These experiments and the theory 19 suggest that 
the charge density in graphene covered by a metal is 
pinned, i.e. not affected by the gate voltage. 

The pinning is much weaker if graphene sheet is oxi- 
dized and the metal-graphcnc contact is resistive. In this 
case the second minimum in the gate voltage dependence 



of conductance can emerge which corresponds to the 
charge density minimum in the metal-covered graphene. 
In Rcf. 19 this mechanism is claimed to be responsible 
for the electron-hole asymmetry observed in Refs. 4 and 
20. We argue, however, that the asymmetry seen in 
Refs. 1, 4, and 5 is entirely due to the contact- induced 
potential in the free-standing graphene. We also regard 
the observed saturation of the conductance for large neg- 
ative gate voltages as the signature of the charge den- 
sity pinning in the metal-covered graphene. The same 
mechanism of the electron- hole asymmetry has been put 
forward recently in Ref. 21. 

The adsorption of graphene on metal substrates has 
been studied by means of the density functional theory 
(DFT) in Refs. 13 and 15. These studies suggest that 
the main effect of metal deposition can be modeled by 
a shift of the Fermi level, i.e. by the chemical doping. 
This shift can be of any sign. Its precise value, /tj, is 
comparatively large due to the small density of states of 
intrinsic graphene, and depends on the workfunctions of 
graphene and the metal as well as on the direct chemi- 
cal interaction. The most recent analysis 15 predicts, e.g. 
/X; = -0.51 eV for Al, -0.40 eV for Ag, -0.43 eV for Cu, 
and 0.21 cV for Au. 

In this paper we refer to the metal-covered graphene 
as the lead, and to the free-standing graphene as the 
sample. The scalar step-function potential in the Dirac 
equation provides the simplest model of the sample-lead 
interface. 6 ' 7 Remarkably, the exact analytical solution 
can also be found if the interface is modeled by an 
exponentially decaying potential. 22 These models, how- 
ever, are not self-consistent and cannot describe the large 
electron-hole asymmetry observed in experiments. 

Microscopic DFT-modeling of the sample-lead inter- 
face has been undertaken in Ref. 18 to investigate the 
contact-induced potential in short (L < 14nm) graphene 
samples, and its effect on the two-terminal conductance. 
In this study a gate electrode is absent and the conduc- 
tance is studied as a function of Fermi energy using clec- 
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trostatic potential for undoped graphenc. The reported 
electron- hole asymmetry is, therefore, of a different type 
than the one occurring with pinned charge density in the 
leads. It is also hard to scale the results of Rcf. 18 to 
realistic system sizes (L > 100 nm). 

Below we propose an effective model, which takes into 
account both the charge density pinning in the leads and 
the self-consistent contact potential in the sample. In 
formulating the model we focus on the two-terminal ge- 
ometry and neglect the resistance of the sample-lead in- 
terfaces. We employ the effective Dirac Hamiltonian and 
take into account the large density of states in a metal. 
The gate electrode is incorporated into the model, while 
the system size poses no serious restriction. The pres- 
ence or absence of a substrate is taken into account by 
an appropriate choice of the permittivity constant. We 
calculate the self-consistent potential in graphene using 
two different approximations, namely the self-consistent 
Hartree (SCH) and the Thomas-Fermi (TF) approxima- 
tion. From the Kubo formula we, then, calculate the 
two-terminal conductance and noise as a function of gate 
voltage. 

We find that the different approximations result in a 
qualitatively similar form of the effective potential. The 
potential penetrates deeply into the sample and reveals 
a power- law decay with the exponent varying from —1 to 
—0.5 depending on the electron concentration. This be- 
havior agrees with the TF analysis of Ref. 14. We demon- 
strate that the slow decay of the potential is responsible 
for the asymmetric gate-voltage dependence of conduc- 
tance and noise, increasing in intensity with the lead dop- 
ing \fxi\. For small positive gate voltages the Dirac point 
crosses the chemical potential twice inside the sample so 
that two n-p interfaces are formed. The slow potential 
decay ensures that the transmission through the inter- 
faces is strongly selective with respect to the momentum 
direction. 23 The electron scattering at the n-p interfaces 
is the reason for the enhanced Fabry-Perot oscillations of 
the conductance as a function of electron concentration 
at moderate positive doping. We note that the oscilla- 
tions transform to the pronounced resonances in the rib- 
bon geometry W < L. 24 When describing gate- voltage 
asymmetries we refer to negatively doped leads, //; < 0. 
For positively doped leads the polarity of the effect is 
reversed. One can further argue that Fabry-Perot oscil- 
lations arc easily distorted by disorder. The presence and 
shape of the Fabry-Perot oscillations can, therefore, be 
used as a direct indicator of the sample quality. 



II. CHARGE AND POTENTIAL PROFILES 

The charge carriers in the graphene sheet (defining the 
xy-plane) are described by the Dirac Hamiltonian 



H = H + V(x) 



H = —ihver • V, 



(1) 



where er is a vector of Pauli matrices and V(x) is the 
spatially dependent effective potential. In Eq. (1) we 



neglect any effects arising from the finite band width and 
assume translational invariance in y. 

Doping of the graphene sheet is caused by the contact 
with the metal electrodes and by applying a voltage to 
the global back gate. We model the lead by graphene 
with the chemical potential ui measured from the Dirac 
point as shown in Fig. 1(a). 

We use the externally compensated charge density 
n{x) in graphene as an input parameter for a self- 
consistent calculation. We focus on a two-terminal setup 
for which the metal-covered and free-standing parts of 
the graphene sheet correspond to |x| > L/2 (leads) and 
\x\ < L/2 (sample), respectively. Accordingly, the com- 
pensated charge density has the spatial dependence 



nix) 



n s , \x\ < L/2 
ni, \x\>L/2' 



(2) 



where n a is proportional to the gate voltage, since the 
distance d g to the gate electrode is typically such that 
the quantum corrections to capacitance can be neglected. 
We assume that the charge density in the lead, ni, is 
pinned, 19 i.e. not influenced by the gate- voltage. The 
screening inside the gate electrode can be neglected for 
L < n s d g , with k s the relative permittivity of the gate 
dielectric. Finally, the inter-valley coupling at the metal- 
graphene interface is also disregarded. 

The relation between the charge densities n s and fii 
and their respective chemical potentials — /j, s and — fj,i 
(cf. Fig. 1(a)) is given by 

»»-,i= / dE Po(E) {(i-e E )(i-f E+f , s t ) - e E f E+ti3 l }, 

(3) 

where f E is the Fermi distribution function, 9 E is the 
Heaviside step function, and 



Po(E) 



■nh 2 v 2 



(4) 



is the density of states (DOS) in infinitely extended ballis- 
tic graphene described by Hamiltonian (1) including the 
spin and valley degeneracy. At zero-temperature Eq. (5) 
becomes 



Hs,l\Vs,l\ 

ttH 2 v 2 



(5) 



Screening in the leads is strongly enhanced by the 
metal electrodes deposited on top of the graphene sheet. 
This effect, which ensures the charge-density pinning, 
is taken into account by means of an additional DOS 
Pm 3> Po{pi)y which is assumed to be energy indepen- 
dent. The transport properties do not depend on the 
value of p m under these conditions. 

The density profile n(x) corresponds to a charge- 
neutral setup with disconnected leads. Bringing the sam- 
ple into electric contact with the leads causes the charge 
redistribution and the band bending shown in Fig. 1(b). 



3 



(b) 




-L/2 



FIG. 1. (Color online) Schematic illustration of the doping in 
the left lead (metal-covered graphene) and the sample (bal- 
listic graphene) if they (a) are not and (b) are in electrical 
contact with each other. The indicated excess charge density 
n e (x) in (b) (hashed regions) corresponds to the semiclassi- 
cal (local) approximation. Quantum mechanical (non-local) 
corrections to n e (x) exist in the gray-shaded regions (Friedel 
oscillations) and in regions where the transverse momentum 
is outside the Dirac cone (evanescent modes). 



plane («bg = 1 for suspended graphene and K^ g w 2.0 
for graphene on SiC>2 with the other side exposed to 
air/vacuum). The distance r(x,y) depends on the "ge- 
ometry" of the graphene sheet. In a planar geometry, 
r = yj ' x 2 + y 2 , one has to regularize the interaction term 
by choosing a finite screening length do, which requires 
that the extension of the leads in x is much smaller 
than ao in order to fulfill the charge neutrality condi- 
tion. Alternatively one can let ao — > oo by considering 
a graphene sheet in the form of a cylinder with circum- 
ference W 3> L. This choice yields the same results for 
V(x) with a computational advantage. 

The self-consistent scheme is closed by relating the 
total charge density n(x) to the potential profile V(x). 
The quantum- mechanical expression for n(x) in the self- 
consistent Hartrcc (SCH) approximation reads 



"scH(a;) = / dEf E _ V ( x ) po(E - V(x)) - f E p(x;E). 

(10) 

The local density of states (LDOS), p(x;E) = 
— (4/7r) ImTr G R (r, r; E), is calculated by solving the 
equation for the retarded Green's function, G R , 



(E + it] - H) G R (r, r'; E) = 5{r - r'), 



(11) 



where r\ is a small positive parameter. We note that the 
dependence of the LDOS on the effective potential is non- 
local. Taking advantage of the translational invariance in 
y we introduce the Fourier transform 



G R (v, r'; E) = -L £ exp[iq (y - y')] G R (x, x';E), (12) 



The excess charge density n e (x) in the entire device is 
defined as 



n e (x) = n(x) + n m (x) — n(x), 



(6) 



where n(x) is the total charge density in the graphene 
sheet relative to intrinsic graphene, and 



n m (x) = (V{x) - Hi) p n 



-L/2 



(7) 



is the additional charge density induced in the metal. 

The potential profile V(x) and the total charge den- 
sity n(x) have to be determined self-consistently. In this 
study we ignore the exchange interaction and restrict our- 
selves to the scalar Hartrcc potential 



V(x) 



Hi 



Vn(x) 



where 



W/2 e -r(x-x',v)/aa 

V H (x)=ahvl dx'n e (x')l dy — — — 

-W/2 r(x-x',y) 



(8) 



(9) 



The coupling constant is given by a = ao/^bg, where 
q = e 2 /4ire hv « (c/v) x (1/137) w 2.2 is the fine 
structure constant in ballistic graphene, and Kbg is the 
effective background dielectric constant in the graphene 



where the summation runs over the discrete values of 
conserved transverse momentum, q n = 2irn/W, with n 
integer (periodic boundary conditions in y direction). 

The boundary conditions in x for the Green's function 
G R (x,x';E) in the channel representation are obtained 
from the exact analytical solution for |x| > £ and \x'\ < £, 
where ^ > L/2 is such that 1^(|£|) — Hi- Selecting the 
decaying solution (G R (x,x';E) — > for x — > ±oo) we 
take the limit 77 — >■ 0. Still, since the spectrum is par- 
tially discrete (some modes are confined to the sample), 
it is necessary to keep rj finite in order to maintain com- 
putational stability of the scheme. We consider the limit 
Hv/W <C r\ <C hv/L such that coherence is preserved on 
the scale hv/rj ;§> L. 

The numerical computation of the Green's function 
G R (x, x; E) for |x| < ^ is demanding since a large energy 
range has to be considered. For hi < Hs the lower energy 
bound at which the integrand in Eq. (10) becomes negli- 
gible is roughly ft-2x (/x s — Hi)- To ensure the charge 
conservation one has to keep the parameters W and rj in 
Eq. (10) identical for both p(x;E) and po(E), so that 



Po(E) = 



hvirW 



5> 



IT] 



y/-(E + irj) 2 + {hvqf 



(13) 
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FIG. 2. (Color online) Contact potentials, V(x), (symmetric 
with respect to x = 0) are calculated in the TF approximation 
(solid lines) and SCH approximation (dashed lines, a = 1) 
for different values of the sample doping, p s , and fixed pi — 
-50 hv/aL (ap m hvL = 100, W/L = 200). 



instead of Eq. (4). The expression (13) converges to 
Eq. (4) in the limit W — > oo, 77 — > 0. 

The relation between charge density and effective po- 
tential is local in the TF approximation, which dramat- 
ically improves the computational efficiency. Replacing 
p(x; E) in Eq. (10) by po{E — V(x)), setting the temper- 
ature to zero, and using the ideal DOS (4) leads to the 
scmiclassical expression for n(x), 



n TF (x) 



V(x)\V(x)\ 
nh 2 v 2 



(14) 



which we regard as the charge density in the TF approx- 
imation. 

The set of Eqs. (8), (6), along with Eq. (10) for SCH 
and Eq. (14) for TF approximations, is solved by means 
of an iterative algorithm. It is run until V(x) reaches 
a self-consistency with accuracy of 10~ 3 relative to \pi\. 
For strong interactions (a <~ 1) and large potential steps 
the procedure requires the use of strong damping in each 
iteration to ensure convergence (the weight of a new it- 
eration is roughly of the order of 10~ 3 ). 

Potential profiles V{x) for different values of p s but 
fixed Hi, calculated at zero temperature, are shown 
in Fig. 2. The dimensionless potential aV(x)L/hv in 
the TF approximation depends on the two parameters 
ap s L/hv and apiL/hv. Such scaling is only approxi- 
mate for the SCH potential. We also note that the TF 
potential transforms exactly as V(x) — > —V(x) under 
the global transformation h(x) — > —h(x), which is not 
the case for the SCH potential. 

Our results apply to large system sizes which are not 
accessible by the DFT models. For instance, for Al 
contacts (pi = — 0.51eV) and L = 250 nm one finds 
piL/hv w —200. The sample width is of minor influence 
as far as W > L. Although the parameter p m determines 
the potential decay in the leads (but not the charge pro- 
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FIG. 3. (Color online) The dotted curve (right axis) shows 
the self-consistent TF potential, V(x), for p„ = and pi = 
-50 hv/L (a = 1, p m hvL = 100, W/L = 200). The corre- 
sponding TF excess charge density n e (x) is plotted with a 
solid line. The dash-dotted curve is the quantum correction 
to the charge density, obtained from the SCH method. This 
correction is split into contributions from energies close and 
further away from the Dirac point as shown by dashed lines. 



file), it has no significant effect on transport properties 
as long as \p s \ < \pi\. 

The validity of the TF approximation is governed by 
the semi-classical criterion \d\(x)/dx\/2TT <C 1, where 
X(x) is the de Broglie wavelength. This is equivalent 
to the condition \dV{x)/dx\ <C V 2 (x) / hv , 14 which is vi- 
olated for the upper curve in Fig. 2 (ap s L/hv = 20). 
The notable difference between the TF and SCH poten- 
tial in this case indicates the importance of non-local 
quantum effects: evanescent modes and Friedel oscilla- 
tions. Evanescent modes strongly increase the charge 
density in a vicinity of the sample-lead interface for 
energies close to the Dirac point. 25 This effect, how- 
ever, is partially compensated in the considered geom- 
etry by the Friedel oscillations. The latter suppress the 
LDOS near the interface for energies far from the Dirac 
point. To illustrate the compensation we plot in Fig. 3 
the quantum correction nscH(^) — ^tf(^) for energies 
\E — V(x)\ < 10 hv/L (dominated by the evanescent 
modes) and \E — V(x)\ > 10 hv/L (dominated by the 
Friedel oscillations). The partial compensation of the 
non-local quantum corrections makes the TF approxi- 
mation reliable even outside its applicability range. Still 
the mentioned discrepancy between the TF and SCH po- 
tentials shows the limitations of local density approxima- 
tions such as TF in positively doped samples. 26 

To make a direct comparison of our self-consistent cal- 
culation to the TF analysis of Ref. 14, we calculate the 
TF potential for a single sample-lead interface at x = 
(fi(x) = hi for x < and h{x) = h s for x > 0). We find 
that the TF potential decays for x > as x~ p . The ex- 
ponent p is given by p = 1/2 for \ V{x)/ pi \ <IC 1 and p = 1 
otherwise (which means that p can vary with position). 
A similar behavior is found in Ref. 14 when the doping in- 
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FIG. 4. (Color online) Potential profiles V(x) for fii = 
— 0.6eV, jU s = 0, a — 2.2, and three different sample lengths 
L — 3.4, 6.8, and 13.6 nm chosen as in Ref. 18. 



side the sample, fi s , is due to the charged impurities. Our 
results for gated graphene arc different since no charge 
density pinning in the leads is assumed in Ref. 14. 

To compare our results to the full DFT treatment of 
Ref. 18 we plot in Fig. 4 the self-consistent TF potentials 
calculated from our model for the same set of parameters 
fii = — 0.6eV, n s = 0, a = 2.2, and L = 3.4, 6.8, and 
13.6 nm. Despite the simplicity of our model (notably, 
the absence of exchange interactions) the potentials agree 
well with those calculated in Ref. 18. 



III. TRANSPORT PROPERTIES 

For a given potential profile one may calculate the 
zero-temperature conductance from the Landauer for- 
mula G = (4e 2 /h)J2 q T q , where the sum extends over 
the transverse momenta q = 2im/W. The shot noise is 
quantified by the Fano factor F = J2 q T q (l - T q ) /J2 q T q- 

The transmission probability T q = T q (E = 0) for 
a given channel at the Fermi energy is related to the 
Green's function by the Kubo formula 27 



T q (E) = (hv) 2 Tr [* x G?(S ) -Z;E)v x G?{Z > -Z;E)] , 

(15) 

where £ is an optional cross-section. We choose £ such 
that V{x) ~ m for \x\ > £. 

The conductance in the TF approximation is shown in 
Fig. 5 as a function of the charge density in the sample, 
n s , for different values of /i/ (using the SCH approxima- 
tion gives very similar results). We note that the conduc- 
tance scales approximately with the parameter fiiL/ahv 
for \ni\L/ahv 3> 1. A similar plot of the Fano factor is 
given in Fig. 6. 

The most evident consequence of the charge transfer 
between the sample and the leads is an electron-hole 
asymmetry in the dependence of the conductance and 
noise on the charge density n s (or the gate voltage). To 
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FIG. 5. (Color online) Zero-temperature conductance as a 
function of external doping n s . The TF potential profiles are 
used. The thin solid lines marked with "a — > 00" correspond 
to step-like potentials, V(x) = Hi for |ar| > L/2 and V(x) = /U s 
for \x\ < L/2, with fiiL/hv given by -50 and -400 for the red 
(light) and blue (dark) line, respectively. 
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FIG. 6. (Color online) Zero-temperature Fano factor as a 
function of external doping n s . As in Fig. 5, the TF potential 
profiles are used (for the thick lines) and the thin solid lines 
correspond to the corresponding step-like potentials. The hor- 
izontal dotted line indicates the asymptotics F = 1/8 for a 
step-like potential with /x; — > ±00. 6 



illustrate the key role of the slow potential decay we plot 
in Figs. 5,6, for comparison, the results obtained from 
the step-function model (V(x) = fi s for \x\ < L/2 and 
V(x) = Hi for \x\ > L/2) with thin solid lines. The 
step- function model also leads to an asymmetry, which is, 
however, negligible in the experimentally relevant regime 
\Hl\L/fw > 1, |/i s | < \m\. 
We find that the conductance is enhanced for neg- 
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FIG. 7. Zero-temperature transmission probability T q as a 
function of transverse momentum q and doping n s (a — 1, 
HiL/hv = -200). 
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FIG. 8. (Color online) Same conductance data as in Fig. 5 
plotted with fixed fj,i and varying L, which is indicated in 
multiples of the sample length L — La for the solid curve. 
Note that the horizontal axis is inverted. 



ative doping since |V(a;)| > \fx s \ in the entire sample 
and suppressed for sufficiently large positive doping since 
|V(x)| < |/i s | in this case. The minimal conductance is 
gradually increasing for large fiiL, while the position of 
the conductivity minimum is shifted towards the positive 
doping. 

Another consequence of the slow potential decay is the 
specific form of the Fabry-Perot oscillations, which are 
only visible on one side of the conductance minimum (in 



the Fano factor one can also observe weak oscillations for 
negative doping). As demonstrated by thin solid lines 
in Fig. 6, a small asymmetry already exists in the step- 
function model with finite doping of the leads. However, 
the oscillation amplitude is strongly enhanced (reduced) 
for positive (negative) doping as a consequence of the 
slow potential decay. The enhancement is strongest when 
the position of the Dirac point, V(x), coincides with the 
chemical potential in a spatially extended region, which 
(as illustrated in Fig. 7) leads to the selective transmis- 
sion of the carries with small transverse momenta at the 
n-p interfaces. 23 

In a particular set of experiments it is the sample 
length, L, rather than the lead doping, that is varied. In 
Fig. 8 we rearrange the data to show the dependence of 
conductance on the sample doping for different sample 
lengths. As expected the contact-induced electron-hole 
asymmetry is more pronounced for shorter samples. The 
shift of the conductance minimum qualitatively agrees 
with recent measurements reported in Ref. 21. 

IV. CONCLUSIONS 

In conclusion we have modeled two-terminal electron 
transport through ballistic graphene samples by taking 
into account the effects of charge transfer at the metal- 
contact/graphene interfaces. Our analysis explains the 
electron-hole asymmetry in the conductance and the 
Fabry-Perot conductance oscillations at positive doping, 
which have been observed in many experiments. 1,2,4,5 
These phenomena are most clearly resolved in Ref. 5 
(Fig. 2 b,c) for two different sample lengths. The period 
8fi s of the oscillations, found by transforming the posi- 
tions of the peaks into fj, s , corresponds to 8[i s L/hv re ir, 
which unambiguously confirms that the oscillations are of 
the Fabry-Perot type. The amplitude of the oscillations 
in this experiment is somewhat stronger than in Fig. 5, 
which is likely due to a modification of graphene by the 
metal leads that has resulted in an increased reflectivity 
of the sample-lead interfaces. 
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